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ABSTRACT 

The evolution of s^arF of about 50 and 12 solar masses has 
been studied through the oxygen-burning phase. For the star of 
50 solar masses it is found that an extensive surface convection 
zone invades deep into the stellar core, and the material in the 
core is mixed into the envelope. When the neutrino loss due to 
the universal Fermi interaction is not taken into account, the re- 
sults suggest that the core mass of a pre -supernova will be close 
to the Chandrasekhar limit. On the other hand, when neutrino loss 
is included, the mixing after the carbon -burning phase is found to 
be negligible, as the time scale of evolution is short compared 
with the time scale of heat transfer even near the bottom of the 
hydrogen -rich envelope. The evolution of a star of about 6 solar 
masses has been studied also. If there is neutrino loss, the 
mixing is found to be strong. Eventually, a pure helium envelope 
develops because of the co-existence of a deepening surface con- 
vection zone and hydrogen- shell burning. 
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I. INTRODUCTION 

The study of stellar evolution in and after the carbon-burning 
phase has been usually done neglecting the hydrogen-rich envelope. 
The present status of such studies is summarized in a recent paper 
by Stothers and Chin (1969). Such studies approximate the evolu- 
tion of stellar cores with a given core mass. However, Stothers 
and Chin (1969) have pointed out that extensive mixing between 
core and envelope occurs in the very advanced phases of evolution 
because of the extensive surface convection zone. They have 
studied very general models by assuming constant radiative flux 
in zones in radiative eq.uilibrium, i.e., mathematically, they 
have solved the stellar structure only as a boundary -value problem. 
Such an assumption is valid in cases when the time scale of heat 
transfer is sufficiently short compared with that of evolution. 

When neutrino loss due to the universal Fermi interaction is 
taken into account, however, the time scale of evolution after 
carbon burning is shorter than that of heat transfer, even near 
the bottom of the hydrogen-rich envelope. In order for a mass 
element in the core to be mixed into the surface convection zone, 
the value of its entropy should be increased up to the value in 
the surface convection zone. Thus, it is important to include 
the equations of energy conservation and of heat transfer. In the 
present study, we have solved the stellar structure and evolution 
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as a combined initial -boundary value problem by means of a 
Henyey-type computational method (Sugimoto 1969 a). 

In the present Note we give only the essentia] assumptions 
and results in so far as they are related to the mixing between 
core and envelope. Other aspects of the evolution, and further 
discussion will be given in separate papers. 

II. HYDROGEN-RICH ENVELOPE AND SURFACE CONVECTION ZONE 

In order to simplify the treatment, the hydrogen-rich envelope 
is replaced by boundary conditions at its innermost layer. In a 
recent study Arnett (1969) has imposed boundary conditions at a 
certain point in the carbon-oxygen core to take into account the 
outer parts of the stars. However, we can formulate the boundary 
cone ■ tions in better approximation at the bottom of the hydrogen-rich 
envelope, because the core interior to it is of a strongly centrally- 
condensed type (Hayashi, H5shl, and Sugimoto 1962j to be referred to 
as ”HHS"). 

We denote values at the bottom of the hydrogen-rich envelope 
with the subscript H^ adding e to indicate the exterior side of the 
chemical composition discontinuity. Other symbols have their usual 
meanings (see, e.g., HHS) For Ug = (4 ■nr^p/M^)^ « 1, the 
following boundary conditions very well take place of an envelope 
solution (HHS): 


k 


G M^p 




^C^)h, 6 ' ^ H, e’ 


(la) 


16 tt a c G .. 

+ 1) ^ u =( -Ju 

rad, H, e V. 5 p k L > H, e 


, k (radiative), 

1 (n + 1) (convective). 


(lb) 


ad, H, e 


where (n + denotes adiabatic polytropic index. In equation (lb), 

the convective (radiative) boundary condition is used when the 
surface convection zone does (does not) reach the bottom of the 
hydrogen -rich envelope. 

To find the depth of the surface convection zone, we formulate 
the so-called "surface condition" of the star in a different way. 
Following equation (4E26) in HHS, the surface condition can be 
approximated as 


L = Lo (M/N^) (T^^^/10^*^ °K)"^^ (e/20)"°*-^. 


( 2 ) 
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HHS gave the values of log lo/l^ as 3*63 for 15*6 % (HHS; ?D1), 
and 3.22 for 4 (HHS; 7D5)* la an adiabatic envelope the 
following relation holds between pressure and temperature, 

P = K (l + ^ “ • " ) exp ^ 1 . (3) 

6 B 

where B denotes the ratio of gas pressure to total pressure, 
including non -degenerate gas and radiation. The integration 
constant K is related to the non-dimensional constant E by 

E = (M<H/k)^'^^K . (4) 

As seen in equation (3), the quantity, - In [(M'H/k)^^^K] , 
is entropy per particle in units of the Boltzmann constant k 
(denoted as ixcr) plus a constant. For the sake of convenience in 

computing gravitational energy release (Sugimoto 1969t), the ionic \ 

! 

entropy per ion is defined as k In + const . , where | 

P^ denotes ion pressure and the const . is selected so that the value , 

of the ionic entropy is the same as that of electronic entropy per 
electron (with quantum statistical zero point) for non -degenerate 
gas consisting of ^°Ne. i 

Eliminating R and from equations (2), and (4), and the | 
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relation, L = n a c R the value of entropy in the convection 

zone is approximately expressed as 

= 8 - 5-1 ; ( 5 ) 

In equation ( 5 ) the additional term involving InCM^'^^E 
•f constant has been replaced by a mean value for the relevant range 
of masses treated in this Note, since we con n'^t solve for the value 
of E in the framework of the present study. The uncertainty in 
equation ( 5 ) is estimated to be ± 1, which comes from uncertainties 
i- Lo and from the terms involving the total mass and E. This 
uncertainty can be ascribed to an uncertainty in log as small 

as ^ 0.2, or to an uncertainty in the value of the core mass of 
about ± 2 per cent as found after computation of the evolution. 

Here, it is interesting to note that the nature of the mixing 
depends only slightly upon the stellar total mass, which may have 
been reduced as a result of mass loss. 

By assuming a value of the core mass and by using boundary 
conditions (l), we solve for the core structui*e and obtain the 

value of the entropy at the bottom of the hydrogen-rich envelope, 
(|jLa)„ . If this is smaller than P'O’* in equation (5), the surface 
conveccion zone does not reach the bottom of envelope; if larger, a 
part of the core should have been mixed into the convective envelope. 


and we then assume a smaller core mass for the next stage of evolution. 
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III. RESULTS AND DISCUSSIONS 

We assume for the initial stage a pure helium core in the 
helium-burning phase. The products of helium burning were assumed 
to be half carbon and half oxygen in mass. For the initial core mass 
we selected 10 5 and I .5 these are discussed in the 

following subsections. 

a) Core of 10 M- 

This core is formed as a result of hydrogen burning of a star 
of about 30 (Stothers 1963 )* In the phase of contracting carbon- 

oxygen core, which follows the helium burning phase, the surface 
convection zone becomes deeper and deeper. Then, the mixing between 
the core and envelope commences at a stage of the phase of contracting 
carbon-oxygen core. Core masses found just after each burning phase 
are as follow: when there is no neutrino loss, 7*^ ^ after carbon 

burning, 6.9 after neon burning, and 6.2 after oxygen burning. 

On the other hand, when there is neutrino loss, the core mass found 
just after the corbon-burning phase is 8.8 After that only 

0.01 ifc> mixed into the lower envelope until the end of the oxygen- 
burning phase. 

In order for a mass element in the core to be mixed into the 
convective envelope, its entropy must be increased by an amount Ao. 

D 

An order of megnitude upper limit estimate of the mass, 6M, which 
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can be mixed diiring a time At , is expressed as 

6L , At 

^ M < (6) 

(k/p,E) T Aa 

In the above relation denotes the energy absorbed in the outer 

part of the core which is in radiative equilibrium, and it is re- 
lated to tlie time scale of heat transfer. Stothers and Chin (1969) 
have discussed the mixing time, i.e., the time scale for establish- 
ing constancy of entropy in a convective enveloi>e. However, relation 
(6) gives a more stringent restriction in limiting mixing. 

We consider first the case with neutrino loss. The lifetime 
from the onset of mixing to the onset of ciirbon burning is 3*1 x 10^ yr, 
and the lifetime of carbon burning is 3*8 x 10^ yr. Thus the core 
mass is essentially determined at the onset of carbon burning. The 
lifetime from the exhaustion of carbon through the exhaustion of 
oxygen is 3 x 10^ yr. If we adopt Aa ~ 10, T ~ 2 x 10^ °K, and 
~ 10^ relation (6) gives 6M ~ 0.01 in agreement with 
the result of numerical comuutation. Mixing after the oxygen- 
burning phase will be even less. The core mass at the pre-supernova 
stage will be about 8.8 and the supernova explosion will be 
triggered by photodisintegration of iron nuclei; this large core mass 
is one of the main differences between the conclusions of the 
present Note and those of Stothers and Chin (1969). 
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On the other hand, when there is no neutrino loss, the life- 
time of the carbon -burning phase is as long as 9*7 ^ 1^0^ yr* and 
the lifetimes of later phases are of the same order of magnitude. 
TSius, condition (6) provides no restriction at all, and mixing will 
proceed as discussed by Stothers und Chin (1969), i.e., the core 
mass of a pre-supernova will be close to the Jhariirasekhar limit, 
and the products of nuclear burnings will be conveyed outwards to 
the stellar surface. 

b) Core of 3 

The DS.SS fraction of the helium core left after the hydrogen- 
burning Is about a quarter for a star of mass between I5 and 4 

(see, e.g., HHS). Thus the stellar mass corresponding a core of 
3 is about 12 M^. There is no mixing Uirough the oxygen-burning 
phase; this result is in essential agreement with that of Stothers 
and Chin (1969)* After the oxygen -burning phase, the same discussion 
will be applicable as that concerning the core of 10 

Core of 1.5 



The evolution of a helium star of I.5 has been computed with- 
out neutrino loss. It is found that corbon -burning begins in a non- 
degenerate carbon oxygen core. Thus, when there is no neutrino loss, 
the evolution of a cere of I.5 will be similar to that of the core 
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of 5 '^lescribed above. No attempt to fit hydrogen-rich envelope 

by boundary conditions (i) 

was made for the case without neutrinos. 

For the case with neutrino loss wc imposed boundary condi- 
tions (l).The stellar mass corresponding to a core of I .5 H. is 
about 6 After the helium has been exhausted in the central 

region, electrons become degenerate because of cooling by neutrinos 
(see e.g., Weigert I966). The central density increases as the core 
contracts. After the central density has become greater than 
4.5 X 10® g cm”®, the surface convection zone invades deeper and 
deeper into the core. In 1.2 x 10^ yr , the core mass has been 
reduced to I.07 Then two carbon-shell flashes occur successively; 

similar flashes have been found by Murai, Sugimoto, Hoshi, and 
Hayashi (1968), and Kutter and Savedoff • Convection develops 

in the carbon -burning shell, which extends from = 0.55 to 
1.0 M 0 * TI*e carbon in this convective shell is almost consumed 
after the two flashes. 

The entropy distribution through the core becomes greatly 
changed by the flashes, and a hydrogen- shell burning becomes active. 
Simultaneously the surface convection zone becomes deeper and deeper, 
xhus, helium produced by the hydrogen- shell burning is mixed through- 
out the envelope, and subsequently the hydrogen in the envelope is 
uniformly consumed. The core is cooled by neutrino loss and the cere 
becomes smaller and smaller in mass. In 7 x 10® yr hydrogen in the 
envelope is exhausted, i.e., the stellar envelope is composed of 
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helium. It is to he noticed that the above lifetime of hydrogen- 
burning in the envelope is comparable to the lifetime of helium 
burning viz., 3*7 x 10® yr. The core mass at this stage is some- 
what uncertain because of ui. certainties in equation (5) for such a 
star; however, its upper limit is estimated to be 0.8 1^. Numerical 
computation has been quitted at this stage of evolution. 

Subsequent evolution of the core will be similar to that of 
Rose (1969), since the difference in mass of the helium envelope is 
not important for the evolution of such a core. Eventually carbon 
will begin to burn at the center when the mass of the carbon-oxygen 
core grows close to the Chandrasekhar limit by the helium- she 11 
burning as studied by Rose (1969) and Faczynski (1969), and carbon 
detonation will disrupt the star as studied by Arnett (1968, 1969)* 

Even if an insufficient amount of carbon has been produced 
as a result of helium burning to result in a carbon -she 11 flash, 
hydrogen burning will eventually begin when the core mass is reduced 
by mixing and the helium envelope will develop. Products of the 
preceding helium burning as well as products of carbon-shell flashes 
are processed in the hydrogen -burning shell at temperature of about 
1.0 X 10® °K and are conveyed outwards into the envelope. In this 
connection it is interesting to point out that hydrogen -deficient 
carbon star R CrB lies near the Cepheid strip in the H-R diagram 
(Searle I961), where a star with a helium surface convection zone 
is expected to lie (Nariai 1969)» 
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